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Abstract-In recent years, morphometric data for Venus and several outer planet satellites have been collected, 
so we now have observational data of complex craters formed in a large range of target properties. We pre- 
sent general inversion techniques that can utilize the morphometric data to quantitatively test various models 
of complex crater formation. The morphometric data we use in this paper are depth of a complex crater, the 
diameter at which the depth-diameter ratio changes, and onset diameters for central peaks, terraces, and peak 
rings. We tested the roles of impactor velocities and hydrostatic pressure vs. crustal strength, and we tested 
the specific models of acoustic fluidization (Melosh, 1982) and nonproportional growth (Schultz, 1988). 
Neither the acoustic fluidization model nor the nonproportional growth in their published formulations are able 
to successfully reproduce the data. No dependence on impactor velocity is evident from our inversions. Most 
of the morphometric data is consistent with a linear dependence on the ratio of crustal strength to hydrostatic 
pressure on a planet, or the factor c/pg. 


INTRODUCTION 

Crater morphometry, the quantitative description of the shape of 
impact craters, has always played a key role in understanding the cra- 
tering process. One of the key arguments used to support the impact 
origin of lunar craters was that they were morphometrically similar 
to terrestrial explosion craters (Baldwin, 1949). Complex impact cra- 
ters, craters with such features as a flat floor, a central peak, and wall 
terraces, have never been created in common geologic materials in the 
lab or with large explosions. At present, only the morphometry of im- 
pact craters on the solid bodies of the solar system can provide data 
on how various target and impactor properties affect complex crater 
formation. Until recently, morphometric data for fresh complex cra- 
ters existed only for the Moon, Mars, and Mercury (e g., Hale and 
Head, 1979, 1980, 1981; Malin and Dzurisin, 1978; Pike, 1977, 
1980a, 1988; Pike and Spudis, 1987; Smith and Hartnell, 1978; 
Wood, 1980; Wood and Andersson, 1978). Unfortunately, these 
bodies provide data for a fairly limited range of target properties. As a 
consequence, past attempts at explaining interplanetary differences in 
craters have necessarily been limited to a forward-modeling approach, 
where a model is considered adequate if it fits the data within error 
tolerances. A well-known example of this approach is Pike's (1980a, 
1988) demonstration that a line with slope 1/g, where g is surface 
gravitational acceleration, fits within error bars the simple-to-com- 
plex transition diameter for the Moon, Mercury, Mars, and Earth. 

The recent addition of crater morphometric data for Venus (Her- 
rick and Phillips, 1994; Sharpton, 1994; Alexopoulus and McKinnon, 
1994; Herrick et ai, 1996) and the icy satellites (Schenk, 1989, 1991) 
greatly extends the range of surface gravities and target compositions 
for which data exist. This additional data makes the inverse approach 
a feasible method for determining the factors controlling interplane- 
tary differences in crater morphometry and for evaluating various 
complex-crater formation models. In the inverse approach, a gener- 
al model is presented and the data are used to invert for any unknown 
parameters in the model. Standard inversion techniques provide a 
structured framework for comparing models, incorporating data 
with errors, and determining ranges of acceptable parameters. Our 
purpose in this paper is twofold. We will attempt to demonstrate the 


advantages of the inverse method for making quantitative evalua- 
tions of models of the cratering process. In the process, we will test 
some general concepts regarding the effects of projectile velocity, tar- 
get surface gravity, and crustal strength on complex-crater formation, 
and we will test two previously proposed models. 

INVERSION TECHNIQUES 

To illustrate what types of data are necessary, what kinds of 
models can be tested, and what can be learned about each model, 
here we briefly review some of the key aspects of linear and non- 
linear inversion techniques (excellent, more detailed reviews are in 
Menke, 1989 and Lines and Treitel, 1984). The first step in using 
inversion techniques to evaluate a model for some aspect of the cra- 
tering process is to quantify the model so that it predicts a set of 
measurable quantities, i.e., measurement d, is a function of some set 
of properties such as surface gravity and impactor velocity. Thus, 
there is a set of measured data points, listed as a vector d, that are 
presumed to be dependent on some unknown set of model param- 
eters, m. In a general sense, this functional dependence can be ex- 
pressed as f(d,m) = 0. All of the techniques discussed in this paper 
require that this functional dependence must be expressed in a way 
that separates model from data, so that 

d = g(m) Eq. (1) 

where g is a vector function. It is easiest to work w ith linear vector 
functions, so that d is just m multiplied by some matrix G, or 

d = Gm Eq. (2) 

As we demonstrate in the specific examples later in the paper, in 
some cases nonlinear g can be made linear through a change in vari- 
able, working in log-log space, or some other mathematical tech- 
nique. If the model and data can be put into the linear form, then 
generally the goal from there is to estimate m from d; in other 
words, invert Eq. (2) so that m is isolated. 

We can begin the model evaluation by looking at the matrix G. 
If G is such that every different m produces a unique d and every 
possible d can be produced by some choice of m, then the model we 
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have chosen is exactly determined. For example, exactly one line can 
always be fit through any two points. In this case, G would be square 
and invertible mathematically. For the exactly determined model, it 
is a foregone conclusion that the data can be exactly produced by a 
particular choice of model parameters, and the only means of evalu- 
ating the model is whether those model parameters and the model it- 
self make sense. Consider as an example a model where the onset 
diameter of central peak craters is a planet's gravity multiplied by an 
arbitrary constant that is different for each planet. Obviously, an ar- 
bitrary constant can be chosen for each planet so that the data are 
always exactly fit; but unless there is some physical basis for the ar- 
bitrary constant, its value and the model are essentially meaningless. 
An exactly determined model can be thought of as transferring the 
data values into a different set of units. 

If G is such that many different m can produce the same d, then 
the inversion for m is said to be underdetermined. For example, an 
infinite number of lines can be fit through a single point. For the un- 
derdetermined case, evaluation of the model must be essentially in- 
dependent of the data. For example, if projectile mass, velocity, and 
impact angle are inputs into a model that predicts the volume of melt 
for a particular crater, then creating values for mass, velocity, and 
impact angle for an incoming meteoroid that result in a match to the 
estimated melt volume says nothing about the validity of the model 
being used. In simple problems, the underdetermined case occurs 
when G has more columns than rows, or m has more elements than d. 
In summary, the validity of the model cannot be assessed with inver- 
sion techniques for the underdetermined and exactly determined cases. 

However, if G is such that only a subset of possible d's can be 
produced by all choices of m, then the model is considered overdeter- 
mined. For example, of all possible sets of three points, only certain 
subsets can be fit by the equation of a line. In the overdetermined 
case, the inversion provides some means of testing the model. A 
model can be ruled out as invalid if a set of model parameters can- 
not be found that satisfactorily reproduces the observed set of data. 
The nature of the inversion is in part based on what criteria define a 
satisfactory reproduction of the data. In this paper, we seek solutions 
that minimize the square of the difference between the model results 
and the data. This is known as a least-squares inversion. In a simple 
least-squares inversion, the model parameters that best reproduce 
the data are of the form (Menke, 1 989) 

m esi = [G T G] -1 G T d Eq. (3) 

Equation (3) gives equal value to each data point; however, it may 
be decided that it is more important to fit certain data points than 
others. In this case, the data can be weighted in the inversion so that 

m est = [G T WG]-‘G T Wd Eq. (4) 

where W is a square diagonal matrix of weighting values for each 
data point. In this paper, we will weight the data by the inverse of its 
standard deviation, so that data points known with the smallest amount 
of error are weighted most heavily. 

Equation (4) and many other linear inversion formulae have the 
form m t , v/ = Md, where M is the inversion matrix. If the data are in- 
dependent of each other so that their covariance is a diagonal matrix 
V where each diagonal element is of , the variance of a data point, 
then an estimate of the model parameter variance is 

V m = MVM t Eq. (5) 

An a priori estimate of the model parameters <m> can be given 
weight in the inversion so that the output model parameters are 


reasonably close to this initial estimate. Combining weighting of the a 
priori estimate with the weighted least-squares inversion gives 

m est = <m> + [GTWG + f^J-'G^ld - G<m>] Eq. (6) 

where W m is a square diagonal matrix individually weighting each 
model parameter and e 2 provides relative weighting between the data 
and the a priori estimate. Equation (6) can be used to address the non- 
linear case where the problem can be defined in terms of Eq. (1) but 
not Eq. (2). If the functions in Eq. (1) vary smoothly, then each func- 
tion can be expanded as a Taylor series where we ignore second and 
higher order terms, so that 

di” + Eq. (7) 

dmj 

where is some initial estimate of a model parameter that the Tay- 
lor series is expanded about. In matrix form, this becomes 

d * G <m> + G’(m-<m>) Eq (8) 

where G' is the matrix of partial derivatives. Setting Ad = d - G<m> 
and Am = m - <m>, and utilizing Eq. (6), it is possible to create an 
iterative method for solving for the best possible estimate of m. In 
this method, known as the Marquardt-Levenberg method (Levenberg, 
1944; Marquardt, 1963), a previous estimate of m and its model 
output are used for <m> and G<m> and the matrix of derivatives is 
used in the inversion, so that the next estimate can be calculated as 

Am = [G' T WG' + /?] -1 G ,T WAd Eq. (9) 

The damping factor/? generally reduces the step size to stabilize the 
inversion. Typically iterations are continued until the step size falls 
below a specified value. 

Once the model parameters are estimated, one way to evaluate 
the model is to compare the prediction error, the difference between 
model results and the data, with the standard deviation a of the data. 
Comparison of different models designed to reproduce the same data 
must take into account the fact that the model with more parameters 
has an inherent advantage in reproducing the data. This comparison 
can be performed with the F test (Menke, 1989). The F value is the 
ratio of Xv f° r each model, w'here 

Xv = l IvZef fo'jFiv^N-M Eq. (10) 

/=l 

e is prediction error, N is the number of data points, M is the number 
of model parameters, and v is the number of degrees of freedom. 
Note that xl becomes smaller with lower prediction errors but 
larger with more model parameters. The ratio, 

F ~xlJ%\ Eq. (11) 

can be compared with standard statistical tables to determine wheth- 
er a significant difference exists between the two models. 

HYDROSTATIC PRESSURE VERSUS STRENGTH 

A particularly illustrative example of the inversion technique is 
to test the role of hydrostatic pressure in the complex crater forma- 
tion process. Many conceptual models of complex crater formation 
envision collapse of a parabolic transient cavity when the hydro- 
static pressure at the cavity's base exceeds some measure of rock 
strength (e g., Hartmann, 1972; Melosh, 1977, 1982; Grieve et al , 
1981). If proportional growth is assumed so that the depth/diameter 
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ratio of a transient cavity is invariant in all cases, then the crater diam- 
eter for collapse will depend on the factor c/pg , where p is surface 
density, g is surface gravitational acceleration, and c is some measure 
of crustal strength. For the terrestrial planets, the crustal strength and 
density are thought to be similar, leading to inverse gravity scaling or 
a dependence on the factor \/g. In particular, Pike (1980a, 1988) was 
able to show that a wide variety of collapse features (e g., inflection 
of depth-diameter data, central peak onset, peak ring onset) on the 
Moon, Mercury, Mars, and Earth can be fit with a line of Mg slope. 
Since that work, the equivalent morphological data has been collected 
for Venus and several icy satellites, and this greatly expands the range 
of available data for all three variables in the factor c/pg. This addi- 
tional data also makes it feasible to use the inversion approach to 
test the importance of the factor c/pg in the complex cratering pro- 
cess. A qualitative, forward-modeling approach to this test was taken 
in Herrick and Phillips (1994), and a preliminary inversion was per- 
formed in Herrick et al. (1996). 

To test the importance of the factor c/pg y we assume that all 
interplanetary differences in crater morphometry are directly con- 
trolled by this factor. If the range of values assumed by this factor is 
small enough, then any dependence can be approximated by an ex- 
ponential function 



where O ^ is an observation of type i for planet j\ and A and b are 
constants. In log-log space this becomes 

(log 0)ji = (log A), + b ,( log c) r b t {[ogpg)j Eq. ( 1 3) 

and a series of equations in the form of Eq. ( 1 ) exists. If is held con- 
stant, then the inversion is linear and Eq. (4) is used for the inversion; 
but if b is allowed to vary, then the nonlinear form of Eq. (9) must be 
used. In this paper, we show results from both cases. 

The particular measurements we used are the following: com- 
plex crater depth, the inflection point of the depth-diameter curve, 
diameter of central peak onset, diameter of peak ring onset, and diam- 
eter of terracing onset. Complex crater depths are functions for each 
planet rather than the single value necessary for Eq. (12). Therefore, 
we used the depth of a 30 km diameter crater, a diameter at which 
craters on all the planets considered are central-peak craters. There 
are at least three distinct measurements that could be considered the 
onset diameter: the diameter of the smallest crater with a particular 
feature, the diameter of the largest crater without a feature, and the 
median diameter, where there are as many craters below the median 
diameter with the feature as there are craters above the median 
without the feature. We refer to these measurements as the minimum, 
maximum, and median onset diameters, respectively. Which type of 
measurement should be used is dependent on both the particulars of 
a model and subjective evaluation of the data. For example, one 
could take the position that the median diameter is most representa- 
tive of what would occur on an ideal planet and that the spread about 
this value represents natural scatter. In this case, the median value 
would be most appropriate, but the scatter should be similar on all 
the bodies so any of the three measurement types would be okay. 
From a different point of view, one might consider, for example, the 
minimum onset of central peaks to be the true onset diameter, and 
larger craters without central peaks have had their peaks covered by 
later crater infilling. Obviously other scenarios can be constructed. 


Table 1 shows the values we used in our various analyses along 
with the appropriate references. The surface gravities (m s -2 ) of the 
planets for which data were collected for this study are: Venus, 8.87; 
Mercury, 3.78; Mars, 3.72; Moon, 1.62; Ganymede, 1.43; Titania, 
0.372; Rhea, 0.285; Ariel, 0.251; Dione, 0.224; and Mimas, 0.079. 
We attempted to maximize the chance of comparing apples with ap- 
ples by generally using data for fresh craters occurring in a single, 
hopefully nonlayered terrain type. For example on the Moon, fresh 
highlands craters were used and on Venus, fresh craters on the vol- 
canic plains were used. Also, crater fields, crater chains, highly 
oblique impacts, and other oddities were not included in the data 
used. Because terrestrial morphological data is both scant and must 
be reconstructed with an implicit model from a heavily eroded 
surface, we have chosen not to use any terrestrial data in our study. 
As stated above, an inversion can be weighted by the error in the 
data. Calculating errors for the depth measurements and the in- 
flection point of the depth-diameter curve is a standard statistical 
problem, and in many cases the error values are given in the stated 
references. However, determining error bars for the onset of a fea- 
ture is not as straightforward. To estimate the error, we first assume 
that our crater data set is a typical sampling of a hypothetical data 
set containing an infinite number of craters. Our sample set of mea- 
surements contains N craters, where N is the number of craters with 
a feature and below the median diameter plus those craters without a 
feature and above the median diameter. An infinite number of sam- 
pled data sets of our fictitious infinite crater set should have a nor- 
mal distribution of median values. If we assume that our observed 
data is a typical sampling of the hypothetical crater population of 
infinite size, then the 1 a error bounds for the median value can be 
estimated by taking the difference of the two data values sam- 
ples away from the median. Using a similar logic for minimum and 
maximum onset diameters, the la error is defined by taking the 
difference of the observed value's diameter and the diameter for sam- 
ple i away from observed value such that 

1-[^J >0.68 Eq. ( 1 4) 

For example, consider the observed data for Venus regarding the on- 
set of central peaks. Figure 1 shows the data between the minimum- 
sized crater with a central peak and the maximum-sized crater without 
one. One hundred ninety-seven craters fall between the minimum- 
sized crater with a central peak, at 8.6 km diameter, and the maximum- 
sized crater without a central peak, 22.6 km in diameter. There are 
27 craters above 14.8 km diameter without a central peak and 27 
craters below 14.8 km with a central peak, making 14.8 km (1.17 in 
log space) the median diameter and N equal to 54. Counting eight 
samples each way (rounding up from -s/54 ) gives an error for the 
median in log space of 0.5 (log 1 5.7 - log 13.5), or 0.033. With N = 
54, i must equal 2 for Eq. (14) to hold. Counting two samples 
away, the error for the maximum is (log 22.6 - log 20.5) = 0.042 
and for the minimum is (log 9.6 - log 8.6) = 0.048. If a reference in 
Table 1 did not provide error estimates or the data necessary to esti- 
mate errors, for the purposes of the inversion we assumed an error 
in log-log space of 0.1, or -25%. 

We performed a series of linear inversions with b fixed at one 
unless otherwise stated. A target density of 3000 kg m 3 was as- 
sumed for the terrestrial planets and moons and 1000 kg m -3 for the 
icy satellites. For the first inversion, we simply tried to fit all of the 
data: rim-floor depth of a 30 km crater, the inflection point in the 
depth-diameter curve, and all three onset diameter measurements for 
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TABLE l . Data and model results for inversions involving hydrostatic pressure dependence. 


Data expressed as log (km) Model Results 



Data 


Error 


Reference 

A 

B 

C 

D 

E 

F 

nlA 

nlA* 

Depth of 30 km crater, Rim- 















Floor (RF) or Terrain-Floor (TF) 

RF 

TF 

RF 

TF 

RF 

TF 

RF 

RF 

RF 

TF 

RF 

RF 

RF 

RF 

Venus 

0.05 

-0.09 

0.14 

0.18 

Sh 

Se, Sh 

0.1 1 

0.22 

0.21 

0.07 

0.19 

0.21 

0.16 

0.21 

Mercury 

0.28 

0.04 

0.04 

0.05 

P3 

P3 

0.32 

0.30 

0.27 

0.05 

0.31 

0.26 

0.32 

0.31 

Mars 

0.22 

0.11 

0.07 

0.31 

P2 

PD, P2 

-0.08 

-0.05 

0.12 

-0.11 

-0.02 

-0.00 

0.01 

0.10 

Moon 

0.47 

0.30 

0.05 

0.08 

P2 

PI, P2 

0.51 

0.56 

0.51 

0.28 

0.52 

0.58 

0.47 

0.41 

Ganymede 

-0.01 

-0.16 

0.09 

0.66 

S2 

S2 

0.06 

-0.08 


-0.12 

-0.04 

-0.02 

0.12 

0.18 

Rhea 

0.40 


0.20 


SI 


0.43 

0.44 



0.45 

0.45 

0.41 

0.37 

Ariel 

0.28 


0.34 


SI 


0.46 

0.38 



0.49 

0.48 

0.42 

0.36 

Dione 

0.40 

0.32 

0.35 

0.48 

SI 

S1.S2 

0.44 

0.46 


0.24 

0.46 

0.54 

0 41 

0.37 

Mimas 

0.52 


0.45 


SI 


0.38 

0.41 



0.40 

0.50 

0.36 

0.34 

d/D inflection point 















Mercury 

0.67 


0.06 


P3 


0.87 

0.86 

0.70 

0.77 

0.85 

0.69 

0.87 

0.87 

Mars 

0.49 


0.09 


P3 


0.46 

0.52 

0.55 

0.55 

0.52 

0.56 

0.45 

0.46 

Moon 

1.04 


0.08 


P3 


1.06 

1.12 

0.94 

1.10 

1.07 

111 

1.06 

1.06 

Ganymede 

0.68 


0.10 


S2 


0.60 

0.48 


0.53 

0.51 

0.58 

0.60 

0.61 

Rhea 

1.09 


0.11 


SI 


0.98 

1.00 


1.01 

0.99 

1.04 

0.98 

0.98 

Ariel 

1.23 


0.10 


SI 


1.00 

0.94 


1.11 

1.03 

1.15 

1.00 

0.96 

Dione 

1.24 


0.17 


SI 


0.98 

1.02 


1.05 

1.00 

1.13 

0.98 

0.98 

Mimas 

1.20 


0.32 


SI 


0.92 

0.97 


0.99 

0.95 

1.09 

0.92 

0.92 

Onset of central peaks 















Minimum 















Venus 

0.93 


0.05 


HP 


0.82 

0.93 



0.89 

0.93 

0.83 

0.83 

Mercury 

1.09 


0.03 


P3 


1.04 

1.00 



1.01 

1.03 

1.04 

1.05 

Mars 

0.52 


0.04 


P2 


0.64 

0.66 



0.68 

0.63 

0.62 

0.61 

Moon 

1.26 


0.05 


SH 


1.23 

1.26 



1.22 

1.25 

1.24 

1.25 

Ganymede 

0.60 


0.10 


S2 


0.77 

0.63 



0.66 

0.57 

0.77 

0.78 

Ariel 

1.00 


0.10 


SI 


1.18 

1.09 



1.19 

1.00 

1.18 

1.14 

Dione 

1.26 


0.10 


SI 


1.16 

1.16 



1.16 

1.15 

1.16 

1.16 

Mimas 

1.00 


0.10 


CM 


1.10 

1.11 



1.10 

1.10 

1.10 

1.10 

Median 















Venus 

1.17 


0.03 


HP 


1.06 

1.12 

1.14 

1.13 

1.10 

1.13 

1.06 

1.06 

Mercury 

1.18 


0.06 


P3 


1.28 

1.20 

1.19 

1.10 

1.22 

1.24 

1.28 

1.28 

Mars 

0.78 


0.10 


P2 


0.87 

0.85 

1.05 

0.88 

0.89 

0.83 

0.86 

0.84 

Moon 

1.43 


0.07 


SH 


1.47 

1.45 

1.44 

1.43 

1.43 

1.45 

1.47 

1.48 

Ganymede 

0.70 


0.10 


S2 


1.01 

0.82 


0.86 

0.88 

0.77 

1.01 

1.01 

Rhea 

1.18 


0.15 


CM 


1.39 

1.34 


1.34 

1.36 

1.26 

1.39 

1.39 

Ariel 

1.18 


0.15 


SI 


141 

1.28 


1.44 

1.40 

1.20 

1.41 

1.37 

Dione 

1.30 


0.11 


SI 


1.39 

1.36 


1.38 

1.37 

1.35 

1.39 

1.39 

Mimas 

1.30 


0.10 


CM 


1.33 

1.30 


1.32 

1.31 

1.30 

1.33 

1.33 

Maximum 















Venus 

1.35 


0.04 


HP 


1.18 

1.28 



1.25 

1.28 

1.18 

1.18 

Mercury 

1.32 


0.03 


P3 


1.39 

1.36 



1.37 

1.38 

1.39 

1.39 

Mars 

1.04 


0.10 


P2 


0.99 

1.01 



1.03 

0.97 

0.97 

0.96 

Moon 

1.70 


0.05 


SH 


1.58 

1.61 



1.58 

1.59 

1.59 

1.59 

Ganymede 

0.85 


0.10 


S2 


1.13 

0.98 



1.02 

0.92 

1.12 

1. 12 

Ariel 

1.30 


0.10 


SI 


1.53 

1.44 



1.55 

1.35 

1.52 

1.49 

Dione 

1.40 


0.10 


SI 


1.51 

1.52 



1.52 

1.49 

1.51 

1.51 

Mimas 

1.54 


0.10 


CM 


1.45 

1.46 



1.46 

1.45 

1.45 

1.45 

Peak ring onset 















Minimum 















Venus 

1.65 


0.00 


HP 


1.65 

1.65 



1.65 

1.65 

1.65 

1.65 

Mercury 

1.86 


0.08 


PS 


1.87 

1.73 



1.77 

1.68 

1.87 

1.88 

Mars 

1.65 


0.10 


w 


1.47 

1.39 



1.44 

1.47 

1.44 

1.42 

Moon 

2.13 


0.11 


PS 


2.06 

1.99 



1.99 

2.04 

2.08 

2.09 

Ganymede 

1.34 


0.10 


S2 


1.60 




1.43 


1.60 

1.59 

Titania 

2.13 


0.10 


S2 


1.83 




1.78 


1.85 

1.97 

Median 















Venus 

1.82 


0.08 


HP 


1.85 

1.96 

1.95 

1.99 

1.92 

1.94 

1.85 

1.85 

Mercury 

2.04 


0.11 


PS 


2.07 

2.04 

2.00 

1.96 

2.04 

1.97 

2.07 

2.08 

Mars 

2.08 


0.14 


P3 


1.66 

1.69 

1.86 

1.74 

1.71 

1.76 

1.64 

1.62 

Moon 

2.34 


0.10 


PS 


2.26 

2.29 

2.25 

2.28 

2.25 

2.33 

2.27 

2.29 

Ganymede 

1.60 


0.10 


S2 


1.80 




1.69 


1.79 

1.79 
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central peaks, peak rings, and terraces. For the icy satellites, we used 
the onset of central pit craters as data for the onset of peak rings, as 
some authors (Melosh, 1982; Passey and Shoemaker, 1982; Schultz, 
1988) have considered central pits and peak rings to be formed by 
similar processes. The form of Eq. (13) is such that only relative 
strengths for the planets can be solved for, and we arbitrarily fixed 
the lunar value of c to be 10(1 in log space). To force the central 


pit onset data to be overconstrained, we assigned Ariel and Titania the 
same strengths. This inversion has 70 data values and 20 parameters 
we are inverting for (/fs and c's in Eq. (13)). Tables 1 and 2 summa- 
rize the results for this inversion, which we will call Model A. Some 
of the fits of model to data are shown graphically in Fig. 2. This in- 
version was surprisingly successful. With the model parameters that 
were inverted for, the forward model is able to fit 51 of the 70 data 


TABLE 1. Continued. 


Data expressed as log (km) Model Results 



Data 


Error 


Reference 

A 

B 

C 

D 

E 

F 

nIA 

nIA’ 

Maximum 















Venus 

1.88 


0.03 


HP 


1.87 




1.88 


1.87 

1.87 

Ganymede 

1.72 


0.10 


S2 


1.82 




1.66 


1.81 

1.81 

Terracing onset: D or d at onset 

D 

d 

D 

d 

D 

d 

D 




d 


D 

D 

Minimum 















Venus 

0.93 

-0.12 

0.02 

0.14 

HP 

HP, Sh 

1.04 




0.09 


1.04 

1.04 

Mercury 

1.26 

0.22 

0.00 

0.04 

SH 

SH, P3 

1.26 




0.21 


1.26 

1.26 

Mars 

0.83 

-0.04 

0.06 

0.07 

P2 

P2 

0.85 




-0.12 


0.83 

0.81 

Moon 

1.20 

0.38 

0.19 

0.05 

SH 

SH, P2 

1.44 




0.42 


1.46 

1.46 

Ganymede 

Median 

1.30 

-0.09 

0.10 

0.09 

S2 

S2 

0.99 




-0.14 


0.98 

0.98 

Venus 

1.11 

-0.06 

0.03 

0.14 

HP 

HP, Sh 

1.13 




0.16 


1.13 

1.13 

Mercury 

1.41 

0.28 

0.10 

0.04 

SH 

SH, P3 

1.35 




0.28 


1.34 

1.34 

Mars 

0.90 

-0.02 

0.10 

0.07 

P2 

P2 

0.94 




-0.05 


0.92 

0.90 

Moon 

1.49 

0.46 

0.03 

0.05 

SH 

SH, P2 

1.54 




0.49 


1.55 

1.55 

Ganymede 

Maximum 

1.60 

0.05 

0.10 

0.09 

S2 

S2 

1.08 




-0.07 


1.07 

1.07 

Venus 

1.29 

-0.01 

0.04 

0.14 

HP 

HP, Sh 

1.36 




0.24 


1.36 

1 36 

Mercury 

1.69 

0.40 

0.23 

0.04 

SH 

SH, P3 

1.57 




0.36 


1.58 

1.58 

Mars 

1.04 

0.04 

0.10 

0.07 

P2 

P2 

1.17 




0.03 


1.15 

1.14 

Moon 

1.62 

0.50 

0.21 

0.05 

SH 

SH, P2 

1.76 




0.57 


1.78 

1.79 

Ganymede 

1.78 

0.12 

0.10 

0.09 

S2 

S2 

1.31 




0.01 


1.30 

1.30 

rms error: data 







0.13 

0.14 

0.08 

0.20 

0.13 

0.14 

0.13 

0.13 

rms error: model 







0.17 

0.13 

0.11 

0.14 

0.14 

0.10 

0.17 

0.17 


References: PI = Pike, 1977; P2 = Pike, 1980a; P3 = Pike, 1988; SI = Schenk, 1989; S2 = Schenk, 1991; Sh - Sharpton, 1994; Se = Sharpton et 
al , 1994; PD = Pike and Davis, 1984; HP = Herrick and Phillips, 1994; CM = Chapman and McKinnon, 1986; W = Wood, 1980; PS = Pike and 
Spudis, 1987; SH = Smith and Hartnell, 1978. 
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Fig. 1 . Raw data for the onset of central peaks for craters on Venus. Graph shows the diameters of craters with and without a central peak. The minimum 
onset diameter is the diameter of the smallest crater with a central peak, or 8.6 km. The maximum onset diameter is the diameter of the largest crater without 
a central peak, or 22.6 km. The median onset diameter is where the number of craters below the median and with a central peak equals the number above the 
median and without a central peak, or 14.8 km. 
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TABLE 2. Model parameters and their formal errors for inversions involving hydrostatic pressure. 


Model parameters 

A 

a 

B 

Model value followed by standard deviation (in log space) 
a C a Da E a F 

a 

nIA 

nIA' 

Log(/f) 















Rim-floor depth, D = 30 km 

3.20 

0.03 

3.24 

0.03 

3.20 

0.04 



3.21 

0.03 

3.21 

0.39 

2.28 

0.47 

Terrain-floor depth, D = 30 km 







2.16 

0.05 







d/D inflection point 

3.74 

0.04 

3.80 

0.04 

3.63 

0.05 

3.79 

0.05 

3.75 

0.04 

4.44 

0.39 

3.50 

1.18 

Central peak onset minimum 

3.92 

0.03 

3.95 

0.03 





3.91 

0.02 

3.12 

0.35 

3.68 

1.37 

Median 

4.15 

0.03 

4.14 

0.03 

4.13 

0.05 

4.12 

0.05 

4.12 

0.03 

3.32 

0.35 

3.92 

1.60 

Maximum 

4.27 

0.03 

4.30 

0.03 





4.27 

0.02 

3.47 

0.35 

4.03 

1.72 

Peak ring onset minimum 

4.75 

0.02 

4.67 

0.03 





4.67 

0.03 

4.97 

0.40 

4.59 

2.22 

Median 

4.95 

0.05 

4.98 

0.05 

4.94 

0.06 

4.97 

0.06 

4.94 

0.05 

5.26 

0.40 

4.79 

2.42 

Maximum 

4.96 

0.04 







4.90 

0.04 



4.81 

2.43 

Terracing onset minimum 

4.13 

0.02 










3.95 

1.59 


Median 

4.22 

0.03 











4.04 

1.67 

Maximum 

4,45 

0.04 











4.27 

1.91 

Terracing onset minimum depth 








3.11 

0.03 





Median 









3.18 

0.03 





Maximum 









3.26 

0.03 





b or a (For models A - F, b is fixed and equal to 1 .0 unless otherwise noted) 















b 




a 


b 

b 

Depth of 30 km crater 







0.7 (fixed) 



0.04 

0.29 

0.67 

0.02 

d/D inflection point 











-0.48 

0.29 

0.91 

0.04 

Central peak onset 











0.61 

0.26 

0.91 

0.05 

Peak ring onset 











-0.18 

0.29 

0.94 

0.05 

Terracing onset diameter 













0.93 

0.05 

Log(c) 















Venus 

1.33 

0.02 

1.40 

0.03 

1.43 

0.05 

1.44 

0.06 

1.40 

0.03 

1.37 

0.04 

1.28 

-7.32 

Mercury 

1.18 

0.02 

1.11 

0.03 

1.12 

0.04 

1.04 

0.06 

1.15 

0.02 

1.04 

0.06 

1.15 

-3.00 

Mars 

0.77 

0.03 

0.76 

0.04 

0.97 

0.06 

0.81 

0.08 

0.81 

0.03 

0.78 

0.04 

0.68 

-12.55 

Moon 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

0.00 

1.00 

1.00 

Ganymede 

0.01 

0.04 

-0.16 

0.05 



-0.10 

0.08 

-0.09 

0.03 

-0.13 

0.06 

-0.05 

-9.81 

Rhea 

-0.31 

0.09 

-0.35 

0.09 



-0.32 

0.10 

-0.31 

0.09 

-0.36 

0.08 

-0.32 

-2.09 

Ariel/Titan i a 

-0.34 

0.05 

-0.46 

0.06 



-0.28 

0.10 

-0.32 

0.05 

-0.37 

0.09 

-0.36 

-2.58 

Dione 

-0.41 

0.06 

-0.43 

0.06 



-0.39 

0.10 

-0.40 

0.06 

-0.37 

0.06 

-0.43 

-2.25 

Mimas 

-0.92 

0.06 

-0.94 

0.06 



-0.90 

0.1 1 

-0.91 

0.06 

-0.87 

0.06 

-0.95 

-4.01 


points within 2 a of the observed value, and 61 of 70 within 3a. The 
rms error (in log space) for the model output is 0.17 vs. 0.13 for the 
actual data, so overall the inversion did not produce an acceptable fit 
to the data. The worst fits of model to data occur for the onset of ter- 
racing and the onset of central pits (treated as peak rings) on the icy 
satellites. The majority of the terracing data is not fit within la, and 
four of the fifteen data values are not fit within 3a of the data. Of the 
four data values for the onset of central pits on Ganymede and Titan ia, 
only one was fit within acceptable error bars. These poor fits suggest 
that the model we used, Eq. (13) with b fixed at 1, was inappropriate 
for these two data sets. It is encouraging that the model outputs for c, 
the strength parameter, seem reasonable. The icy satellites as a group 
have strengths an order of magnitude less than those for the terrestrial 
planets. Mars, which may have water in the crust, has a relative 
strength half that of the Moon. Mercury, which superficially appears 
to most resemble the lunar highlands, has a value of c that is 50% lar- 
ger than the Moon's. Venus has a dry crust that may not be as highly 
fractured as the lunar or Mercurian crust, and its calculated relative 
strength is roughly double the Moon’s. For the icy satellites, Ariel, 
Rhea, Dione and Titania have strengths roughly half that of Gany- 
mede’s, and the strength of Mimas is about one-tenth that of Gany- 
mede's. 

There are many logical variations on this first model that can be 
tested, and we summarize a few here. We ran another inversion iden- 
tical to Model A, except it did not include the terracing data or the 


central pit onset data. The results for this inversion. Model B in 
Tables 1 and 2, were similar to those for Model A and indicate that 
including the terracing and central pit data did not adversely affect 
the rest of the inversion in Model A. The Model B results have a 
slightly lower rms error than the data because the terracing and cen- 
tral pit data were excluded, and all but three of the data points were 
fit within 3a. For the test labeled Model C in Tables 1 and 2, we used 
only data for the rocky planets and the following measurements: 
rim-floor depth of a 30 km crater, the inflection point in the depth- 
diameter curve, and the median onset diameters for central peaks and 
peak rings. The input for Model C was a small subset of that for Mod- 
el A, but the output results were remarkably similar, and the same 
general trend in planetary strengths was found. Without the inclu- 
sion of icy satellites, the rms error for the data was substantially lower 
than the error for the data used in Model A. Thus, even though Mod- 
el C had a lower rms model output error than Models A or B, the 
model output does not have an rms error lower than the observed data. 
All of the data were fit within 3a with Model C. 

Recently McKinnon et ai (1997) stated that a g~ 0J trend was 
more logical for the depth trend for a single diameter, and they also 
preferred to use terrain-floor depth data rather than the rim-floor 
depths that are commonly used. We ran a set of four tests where we 
ran through the possible combinations of using terrain-floor depth vs. 
rim-floor depth and forcing a (c/pg) x vs. ( c/pg ) 0 - 7 trend for the depths. 
The measurements used for the inversions were depth of a 30 km era- 
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FIG, 2. Graphical representation of the data and model results for Model A in Tables 1 and 2. Data is shown as x’s (horizontal lines for median peak ring 
onset) with vertical lines representing 1 a error bars. Model results are shown as diamonds (squares for median peak ring onset). Data and values are plotted 
vs. target surface gravity. Data and results for Mars have been plotted at a slightly lower than actual surface gravity to allow the points to be distinguished 
from the Mercurian data and results. 


ter, the inflection point in the depth-diameter curve, and the median 
onset diameters for central peaks and peak rings. These four inver- 
sions all yielded similar results that were also consistent with Mod- 
els A and B. Figure 3 compares, for all four inversions, the calculated 
model trends vs. the depth data. The results for one of these, the one 
using terrain-floor depths and forcing a (c/pg) 07 depth trend, are 
shown as Model D in Tables 1 and 2. Negligible differences were 
found between model fits using terrain-floor vs. rim-floor depths, 
but the relative strengths of Venus, Mercury, and Mars were slightly 
lower if terrain-floor depths were used. Almost identical model fits 
were found between models using a (c/pg) 07 vs. a (c/pg) 1 trend, and 
the primary model parameter that changed between these two models 
was the constants in Eq. (12) in response to the forced change in b. 

We also performed several nonlinear inversions similar to mod- 
els A— D but allowing b to be a variable for each measurement type. 
For the nonlinear inversions, an initial estimate of model parameters 
must be given as a starting point for the inversion, and damping of 
individual steps may be important. The nonlinear inversion can be 
thought of as starting at a specified point on a surface of model solu- 
tion errors and using local slopes to end up in a local minimum that 
minimizes error. Ideally there is only one local minimum that is the 


absolute minimum, and that minimum is reached from any starting 
point. Unfortunately, as we discuss below, that is not the case for the 
inversions performed here. A judicious choice of the damping factor 
will get the inversion to a local minimum in a few steps without over- 
shooting that minimum. We found that the most effective technique 
was to start with minimal damping and automatically increase damp- 
ing as necessary to ensure that each step produced a model output with 
a lower error than the previous step. We ceased iterating when the av- 
erage step became less than 5 x 1 0 -4 the average model parameter. 

Model nlA in Tables 1 and 2 shows the results for a nonlinear 
inversion that used the results of linear Model A as a starting set of 
parameters. We allowed b to vary, but required b to be the same for 
the minimum, maximum, and median of each measurement type. 
Thus, Models A and nlA are fitting the same data, but Model nlA is 
nonlinear and has five additional parameters (the b's). This particu- 
lar inversion converged to a minimum error solution with fairly minor 
changes in the model parameters from the starting values. In particu- 
lar, the calculated b values stayed near 1.0 except for the 30 km 
crater depths, which inverted for a slope of 0.69. This would seem 
to validate the linear dependence on c/pg except for the (c/pg) 07 de- 
pendence on depth suggested by McKinnon et al. (1997). However, 
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Fig. 3. Crater depth vs. target gravity for 30 km central peak craters on six 
different planetary bodies. Results are from identical inversions (Model D in 
Tables 1 and 2) using the four possible combinations of rim-floor depths or ter- 
rain-floor depths and forcing a g~ 01 or g~ l 0 trend. The top plot shows results 
vj. data for inversions using terrain-floor depths with the two gravity trends 
and the bottom plot shows similar results for inversions using rim-floor depths. 
Data is shown as x’s with vertical lines representing 1 a error bars, and model 
results are shown as diamonds (g -1 0 trend) and squares (g~° J trend). 


there was less than a percent difference in data error between the out- 
puts of the linear and nonlinear models. In other words, allowing h 
to vary essentially gained nothing in terms of the overall model being 
able to fit the data. 

We also found that the solution set of model parameters is non- 
unique. Model nlA* in Tables 1 and 2 shows the results from an in- 
version with a different set of starting parameters. In this case, the 
solution evolved to a set of functions with nearly flat slopes ( b ~ 0) 
and still achieved the same fit of model results to observed data. Two 
characteristics of the observed data account for the nonuniqueness of 
inversion results. First, when plotted on a log-log graph of measure- 
ment value vs. surface gravity, most of the observed data sets are 
roughly parallel. Proportionality constants (the A' s) can be solved for 
that have the effect of collapsing several data sets into one line. Sec- 
ond, we expect the strengths of the planetary bodies to have roughly 
the same trend as the gravity. The rocky planets have higher surface 
strengths and higher surface gravities than the icy satellites. There- 
fore, higher or lower slopes (the b's) can be counteracted by less or 
more dramatic strength variations. Additional nonlinear inversions 
were performed corresponding to the data sets used for some of the 


other linear inversions, and results were similar to those obtained in 
Models nIA and nlA'. 

F tests showed that in no case was the model error reduced by 
using the nonlinear inversions more than could be attributed to simply 
adding parameters to the model. In other words, a simple linear de- 
pendence on the factor c/pg works pretty well at matching most of 
the data we used except for the onset of terracing data. While the 
nonlinear inversion results were not an improvement over the linear 
model, they also showed that a linear dependence is not the only ac- 
ceptable model for fitting the data. 

Model A shows that terracing onset data do not reasonably fit a 
c/pg dependence. Some previous works ( e.g.. Pike, 1980b; Herrick 
and Phillips, 1994) have suggested that the poor fit for terracing on- 
set might be because terracing is a late-stage process that occurs after 
final crater depth has largely been determined. To test this idea in a 
primitive way, we assume that the complex crater depth d, rather than 
diameter, for the onset of terracing is dependent on the factor c/pg. 
We performed another inversion with Model A modified by sub- 
stituting crater depth for crater diameter in the data for onset of ter- 
racing. These results are shown as Model E in Tables 1 and 2. The 
reduction in error is dramatic. For Model E, 58 of 70 model param- 
eters are fit within 2 o and 66 of 70 are fit within 3 o. Overall the 
rms model error is < 1 0% greater than the rms data error. Thus, our 
inversion results support the hypothesis that terracing occurs when 
final crater depth exceeds a constant multiple of c/pg. 

ACOUSTIC FLUIDIZATION 

A more specific model of complex crater formation that has been 
proposed is the acoustic fluidization model of Melosh (1982). The 
basic premise of the acoustic fluidization model is that a hemispheri- 
cal region containing the transient cavity becomes fluidized by 
acoustic noise. The fluidized material behaves as a Bingham plastic, a 
material that behaves as a solid below a yield stress c but flows as a 
fluid when differential stresses exceed this strength. Melosh (1982) 
assumes a fixed viscosity rj for the fluidized material, which then be- 
haves as a damped, oscillating wave that freezes when stresses fall 
below c. The shape of the wave when it freezes determines whether 
the crater has a central peak, a peak ring, or multiple rings. Initial col- 
lapse occurs if 

pgH / c > 5 Eq. (15) 

where H is the depth of the transient cavity, and collapse becomes 
deep-seated if the inequality is >8 to 15. Melosh (1982) assumes a 
depth/diameter ratio of 1 :5 for the transient cavity, so collapse oc- 
curs when 

pgD/c> 25 Eq. (16) 

The break in the depth-diameter function is considered by Melosh 
(1982) to be the best estimate of when collapse occurs, so the diam- 
eters of these breaks and Eq. (16) as an equality can be used to 
invert for c. When collapse is deep-seated, the final depth d of the 
crater should be a constant independent of diameter but determined 
by c/pg so that 

d = Kc/pg Eq. (17) 

where K is a constant that should be between 8 and 15. We can use 
the depth of a 30 km diameter crater and Eq. (17) to invert for K and 
c. Through a series of approximations, Melosh (1982) determines that 
the number of zero crossings that occur in the oscillating wave is, 
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Table 3. Results from acoustic fluidization inversions. 



Approach 1 
data as log (m) 
Data error Results 

Approach 2 
data as In (m) 

Data error Results 

Parameters 

App 1 
(log) 

App. 2 
(In) 

d/D inflection point 







Yield stress (Pa) 



Mercury 

3.67 

0.06 

3.73 

8.45 

0.06 

33.60 

Venus 

6.53 

26.20 

Mars 

3.49 

0.09 

3.62 

8.04 

0.09 

34.10 

Mercury 

6.37 

25.20 

Moon 

4.04 

0.08 

3.97 

9.30 

0.08 

35.00 

Mars 

6.26 

26.00 

Ganymede 

3.68 

0.10 

3.57 

8.48 

0.10 

23.60 

moon 

6.24 

25.70 

Rhea 

4.09 

0.11 

4.05 

9.43 

0.11 

14.90 

Ganymede 

5.33 

15.10 

Ariel 

4.23 

0.10 

4.19 

9.74 

0.10 

11.70 

Rhea 

5.10 

5.45 

Dione 

4.24 

0.17 

4.17 

9.77 

0.17 

14.60 

Ariel 

5.19 

1.96 

Mimas 

4.20 

0.32 

4.13 

9.68 

0.32 

21.20 

Dione 

5.12 

4.87 

Depth of a 30 km crater 







Mimas 

4.63 

1 1 .50 

Venus 

3.05 

0.14 

3.05 

7.01 

0.14 

10.20 

Viscosity (Pa s) 



Mercury 

3.28 

0.04 

3.25 

7.54 

0.04 

8.60 

Venus 

9.34 

1.1*10 4 

Mars 

3.22 

0.07 

3.15 

7.42 

0.07 

9.58 

Mercury 

9.17 

9.5*10 3 

Moon 

3.47 

0.05 

3.49 

7.98 

0.05 

8.75 

Mars 

8.37 

3.2*10 3 

Ganymede 

2.99 

0.09 

3.09 

6.88 

0.09 

-0.78 

moon 

9.56 

1 .4x 10 4 

Rhea 

3.40 

0.20 

3.57 

7.84 

0.20 

-11.40 

Ganymede 

7.86 

8.3* 10 2 

Ariel 

3.28 

0.34 

3.72 

7.55 

0.34 

-14.60 

Rhea 

7.14 

8.6x10' 

Dione 

3.40 

0.35 

3.70 

7.84 

0.35 

-12.00 

Ariel 

7.63 

8.0x10' 

Mimas 

3.52 

0.45 

3.66 

8.11 

0.45 

-5.61 

Dione 

7.28 

1.5xl0 2 

Median central peak onset 


N x N for cent, peak onset 

Mimas 

6.67 

9.3x10' 

Venus 

4.17 

0.03 

4.28 

1.00 

0.15 

7.7 x 10 2 

K 

0.92 

-9.02 

Mercury 

4.18 

0.06 

4.44 

1.00 

0.30 

4^ 

X 

© 




Mars 

3.78 

0.10 

4.45 

1.00 

0.59 

-0.48 




Moon 

4.43 

0.07 

4.$1 

1.00 

0.35 

0.17 




Ganymede 

3.70 

0.10 

4.42 

1.00 

0.62 

-2.68 




Rhea 

4.18 

0.15 

4.18 

1.00 

0.90 

1.11 




Ariel 

4.18 

0.15 

4.18 

1.00 

0.90 

-0.16 




Dione 

4.30 

0.11 

4.30 

1.00 

0.59 

2.09 




Mimas 

4.30 

0.10 

4.30 

1.00 

0.52 

-1.96 




Median peak ring onset 




N x N for peak 

ring onset 




Venus 

4.82 

0.08 

4.23 

4.00 

0.71 

7.50 




Mercury 

5.04 

0.11 

4.23 

4.00 

1.06 

77.40 




Mars 

5.08 

0.14 

3.68 

4.00 

1.52 

51.90 




Moon 

5.34 

0.10 

4.61 

4.00 

0.97 

36.00 




Ganymede 

4.60 

0.10 

3.88 

4.00 

1.00 

50.00 




rms error: data 



0.17 



0.56 




rms error; model 



0.43 



21.83 





See Table 1 for appropriate references for data values. 


N = — 

71 


3 Igp 2 !?' 


V 




In 


2V3 


pgH a 
c L 


Eq. (18) 


where L is the radius of the fluidized region and a is the radius of 
the transient cavity. We use Melosh’s (1982) assumptions that L is 
the also the radius of the final crater and that it is equal to 1.5 a. 
Equation ( 1 8) becomes 


1 

-ii 


define a third constant 

Ntt 

71 

L V 

t c J 

*3 = 

In (k 2 D/c) 


We would like to use our measurements and Eqs. (16), (17), and 
(19) to invert for each planet’s viscosity (in Pa s) and yield stress (in 
Pa). However, Eq. (19) is not in the form required by Eq. (1) with 
data and parameters separated. It is impossible to isolate the data 
values of the onset diameters D in Eq. (19) from the model param- 
eters. We tried two approaches to getting around this problem. In 
the first approach, we held one of the two terms in Eq. (19) fixed dur- 
ing each inversion step and then recalculated that term between steps. 
We chose to hold the second term fixed so that for each step we could 


Eq. (20) 


k\ = 


ki - 


32 

2 pg 

45V3 


Eq. (19) 


For a given planet, N = 1 at the onset diameter of central peaks and 
N = 2 at the onset diameter of peak rings. 


so that now Eq. (19) can be rewritten as 

d 3 = M 2 A 1 

or in log-log space 


log D = - (log lc 3 + 2 log rj - log k , ) 


Eq. (21) 


Eq. (22) 



140 


R. R. Herrick and S. N. Lyons 


In summary, for the first approach we used log-log versions of Eqs. 
(16) and (17) along with Eq. (22) to invert depth-diameter transitions, 
depth data, central peak onsets, and peak ring onsets for the yield 
stress and viscosity parameters. The inversion has 3 1 observations 
to invert for 19 parameters. 

In the inversion, k 3 must be calculated before each step, and this 
value depends on the value of D used. We ran two inversions, one 
where the observed diameters were used and the other where D is cal- 
culated from Eq. (22) using parameters and model results from the 
previous iteration. There were no significant differences in the results 
from either approach, and Table 3 summarizes the results where the 
observed diameters were used to calculate k$. The inversion was un- 
satisfactory in that no set of model parameters was able to success- 
fully match the data within the data error bars. When taken back 
out of log space, the data has an rms error of -50% while the model 
produces errors of— 170%, or the model results are typically a factor 
of 2.7 greater or less than the observed data. The principal failure of 
the model is that no viscosity values could be found so that the 
model could match both central peak onsets and peak ring onsets. Ex- 
amination of Eq. (18) shows that to first order, the model predicts that 
the ratio of onset diameters for peak rings vs. central peaks should 
be 2 2/3 , or -1.6, while observed ratios range from about 4 to 20. 
Consequently, the best the model can do is choose viscosity values 
that produce onset diameters intermediate between the observed val- 
ues for central peaks and peak rings. 

In the second approach, we used Eq. (18) directly with N instead 
of D as the data value we were solving for. We squared Eq. (18) to 
stabilize the inversion. The data values D were treated as constants 
in the inversion, but the percent error in D was used as the percent 
error in A a . Because Eq. (18) involves a natural log term, we used 
ln-ln versions of Eqs. (16) and (17) to invert the depth-diameter 
transitions and depth data. As the results in Table 3 show, this sec- 
ond approach was even less successful than the first approach. In the 
first approach, the yield stresses c were calculated independently of 
r] and were not affected by the problem with the ratio between onset 
diameters for central peaks and peak rings. In the second approach, 
the only way the inversion could produce a large difference in N 
with small differences in D was to choose large strength values that 
made the natural log term in Eq. (19) at or below zero where the 
natural log function is rapidly varying. 

Melosh (1982) obtained a fit of model to data by allowing the 
viscosity to be different for central peak onset and peak ring onset. 
In this case, the inversion is exactly determined so that the viability 
of the model must be evaluated by some other means than the mod- 
el's ability to fit the data. A discussion of the physical viability and 
the geologic evidence for acoustic fluidization in a crater is beyond 
the scope of this paper. It is interesting to note that the basic mathe- 
matics in the acoustic fluidization model should be appropriate for 
any conceptual model of complex crater formation that involves Bing- 
ham fluid collapse of a hemispherical shell of material. Any such 
model must have the property of either decreasing L relative to a with 
increasing crater diameter or increasing rj with increasing diameter. 
This would seem to rule out collapse with the hemispherical cavity 
defined by exposure to a specified shock pressure, as that would pre- 
dict L increasing relative to a (e.g., Melosh, 1989). 

NONPROPORTIONAL GROWTH AND IMPACT 
VELOCITY EFFECTS 

Some authors have advocated that the onset and formation of cer- 
tain properties of complex craters are a result of nonproportional 


growth or a change in the shape of the excavation cavity with in- 
creasing crater diameter. In particular, Schultz (1988) has used ex- 
trapolation of small-scale impact experiments to argue for shallowing 
of transient crater diameter above a crater diameter dependent on 
v 1 17 /g, where v is impactor velocity. He estimates that both the 
depth-diameter transition and the onset of peak rings should be pro- 
portional to the same parameters, so that 


1.17 


D i\,l ~ ^1,2 " 


v 0.22 


Pt_ 

Pp 


\ 0.83 




Jet 

\pt 


0.5 


- 0.83 


Eq. (23) 


where s is sound speed, the subscripts t and p denote target and pro- 
jectile, and the subscripts 1,2 indicate the equations for depth-diam- 
eter transition and peak ring onset, respectively. Schultz (1988) also 
includes an equation for peak ring diameter relative to crater diam- 
eter, but the predicted and observed interplanetary variations are less 
than the error bars on the data, so we did not use that data in our 
inversion. For a given planet, we expect that the highest impact velo- 
cities will generally be associated with cometaiy impactors that we 
expect to have relatively low densities and sound velocities. We 
can carry this logic into the inversion by associating the maximum 
possible impact velocity with the diameter of the largest crater with- 
out a peak ring, the minimum possible impact velocity with the diam- 
eter of the smallest crater with a peak ring, and the median impact 
velocity with the median onset diameter of peak rings and the depth- 
diameter transition. 

For the inversion, we estimated all parameters in Eq. (23) except 
for k\ 2 and inverted for these two constants in log-log space. For the 
Moon, Mercury, and Mars, we assumed a target density of 2900 kg 
m -3 and sound speed of 3 km/s; for Venus, we assumed a target den- 
sity of 2900 kg m -3 and sound speed of 4.5 km/s; and for the icy 
satellites, we assigned a target density of 1000 kg itt 3 and sound 
speed of 1.6 km/s. For the terrestrial planets, we assumed that the 
minimum and median velocity impactors were stony asteroids and 
assigned them a density of 2500 kg nr 3 and sound speed of 3 km/s, 
and the maximum velocity impactors were assumed to be cometary 
and were assigned a density of 1000 kg itt 3 and sound speed of 1 
km/s. For the icy satellites, all impactors were assumed to be come- 
tary. 

Table 4 shows the impactor velocities used in the inversion, the 
constants k\ 2 resulting from the inversion, and the model results com- 
pared to the data. The minimum impactor velocity is simply the es- 
cape velocity of the target body. The maximum impactor velocity' is 
also calculated theoretically using the orbital and escape velocities 
of the sun and the pertinent planets and moons. Median impactor 
velocities for the inner planets are from Schultz (1988); and for 
Ganymede, the Shoemaker and Wolfe (1982) value for short-period 
comets is used. For the remainder of the icy satellites, we use as a 
median velocity the simple approximation of >/3 times the orbital 
velocity of the planet orbited by the satellite. This approximation 
assumes an impacting population of long-period comets and ignores 
the minor acceleration provided by the gravity wells of the planet 
and satellite. All of the velocity values used are speculative, but the 
general trend they define should be accurate. Other than approxi- 
mately defining the general trend by virtue of the inverse gravity 
term in Eq. (23), the model results simply do not match the data. The 
rms model error of 0.65 corresponds to the model results typically 
being a factor of 4.5 greater or less than the data. 
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Table 4. Results from nonproportional growth inversions. 


Density and velocity parameters used in inversion 

Impactor velocities and corresponding impactor density and sound speed Target density 



Min. v 
km/s 

L P/? -3 

kg m 3 

s p 

km/s 

Med. v 
km/s 

. p P_ 3 

kg m 3 

s p 

km/s 

Max. v 
km/s 

. P P 3 
kg m 3 

S P 

km/s 

and sound speed 

Pt . s i 

kg m 3 km/s 

Venus 

10.4 

2500 

3 

27 

2500 

3 

86 

1000 

1 

2900 

4.5 

Mercury 

4.3 

2500 

3 

34 

2500 

3 




2900 

3.0 

Mars 

5.0 

2500 

3 

16 

2500 

3 




2900 

3.0 

moon 

2.5 

2500 

3 

22 

2500 

3 




2900 

3.0 

Ganymede 

2.7 

1000 

1 

16 

1000 

1 

48 

1000 

1 

1000 

1.6 

Rhea 




17 

1000 

I 




1000 

1.6 

Ariel 




12 

1000 

1 




1000 

1.6 

Dione 




17 

1000 

l 




1000 

1.6 

Mimas 




17 

1000 

1 




1000 

1.6 

Titania 

0.8 

1000 

1 







1000 

1.6 


Data and inversion results expressed as log (km) 

Data Error Results Data Error Results 


d/D inflection point 


Mercury 

0.67 

0.06 

0.70 

Mars 

0.49 

0.09 

0.32 

Moon 

1.04 

0.08 

0.85 

Ganymede 

0.68 

0.10 

0.66 

Rhea 

1.09 

0.11 

1.39 

Ariel 

1.23 

0.10 

1.26 

Dione 

1.24 

0.17 

1.49 

Mimas 

1.20 

0.32 

1.94 

Peak ring: minimum onset 
Venus 

1.65 

0.00 

1.65 

Mercury 

1.86 

0.08 

0.82 

Mars 

1.65 

0.10 

0.91 

Moon 

2.13 

0.11 

0.91 

Ganymede 

1.34 

0.10 

1.39 

Titania 

2.13 

0.10 

1.34 


Peak ring: median onset 


Venus 

1.82 

0.08 

2.13 

Mercury 

2.04 

0.11 

1.87 

Mars 

2.08 

0.14 

1.50 

Moon 

2.34 

0.10 

2.02 

Ganymede 

Peak ring: maximum onset 

1.60 

0.10 

2.28 

Venus 

1.88 

0.03 

3.27 

Ganymede 

1.72 

0.10 

2.84 

rms data error 



0.12 

rms model error 



0.65 

Parameters 




log(t|) 



-3.69 

log(* 2 ) 



-2.06 


We also tried a model that tested a more general dependence on 
impact velocity. We modified the linear form of Eq. (12) to allow for 
a possible exponential dependence on velocity, or 

Ojj = — v‘j' Eq. (24) 

PjSj 

where a is an unknown exponential. In log-log space, this expres- 
sion becomes 

(logO) ;V = (log .4); + (logc); -(logpg)y +c / (logv) / Eq. (25) 

which is still linear if we assume we know the impact velocities. We 
ran an inversion with the same data set as Model B but with addi- 
tional parameters defined by the last term in Eq. (25). We used only 
median impact velocities and had a separate a for depth of a 30 km 
crater, the d/D inflection point, onset of central peaks, and onset of 
peak rings. The same a was used for the minimum, median, and maxi- 
mum onset diameters, so a total of four parameters were added to the 
inversion of Model B. 

Our initial inversion was unstable. We decided to perform another 
inversion with some minimal weighting given to an a priori estimate 
to stabilize the solution. For a priori estimates, we used the model 
parameter results of Model B and varied the a values from 0 to 2. 
We found that providing a uniform weighting of the a priori esti- 


mate of -0.1% was all that was necessary to stabilize the solution, 
and the model results using an a priori estimate had an rms error with- 
in 3% of the unconstrained inversion. Model F in Tables 1 and 2 
shows the results with an a priori estimate of 1.0 for the a values. 
In our solutions, only the a value for onset of central peaks consis- 
tently had an absolute value significantly different from zero, and an 
F test comparison between Models F and B indicates that the 
reduction in error from Model B to Model F is not statistically 
significant. We conclude that there is no evidence for an impact velo- 
city dependence for the measurements used in Models B and F, at 
least not in the form of Eq. (25). We emphasize that our results do 
not imply that varying impact velocity has no effect on the volume 
of the excavation cavity but merely support the notion that the shape 
of the excavation cavity is velocity-independent. The observables we 
used in our inversions are for craters formed by hypervelocity im- 
pacts in the gravity-dominated regime, and, thus, our results do not 
apply to differences between subsonic and hypervelocity impacts or 
to craters formed in the strength-dominated regime. 

DISCUSSION 

We have inverted crater morphometry data to test general and spe- 
cific models for complex crater formation. A general dependence on 
hydrostatic pressure v$. strength, or the factor c/pg, works remarkably 
well for most of the data. The trend in planetary crustal strengths that 
is inferred from the inversions is consistent with our expectations: Ve- 
nus has the highest crustal strength, Mercury' and the Moon have 
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0.01 0.10 LOO 10.00 0.01 0.10 LOO 10.00 

Target Surface Gravity Target Surface Gravity 



Target Surface Gravity Target Surface Gravity 


Fig. 4. Similar to Fig. 2 but using model results from Model E in Tables 1 and 2. 


slightly lower strengths because of a more fractured crust, Mars has 
a stilt lower strength perhaps because of water in the crust, and the icy 
satellites are an order of magnitude weaker than the terrestrial planets. 
The acoustic fluidization model of Melosh (1982) did not fit the data 
unless the number of model variables was increased until the model 
was even determined; for even determined models, the fit of model 
to data is guaranteed and is not an indicator of the model's validity. 
The nonproportional growth formulation of Schultz (1988) did not fit 
the data, and no dependence on impact velocity in general was found. 

A dependence on the factor c/pg does not reproduce all of the 
data. The diameter onset of terracing does not follow' a c/pg trend but 
the depth onset of terracing does (Model A vs. Model E). This may 
indicate that terracing is a late-stage process that occurs after the final 
depth of the crater has been determined, an idea first suggested by 
Pike (1980b) using only data from Mars, the Moon, and Mercury. 
While the central pit onset data did not seem to fit the peak ring onset 
data trend in Model A, in Model E the slight change in relative crustal 
strengths results in three of four model results for central pit onset 
falling within error bars for the data. Figure 4 shows graphically the 
model results vs. data for model E. Perhaps what stands out most in 
Fig. 4 is that Mars consistently does not fall on the trend produced 
by the other planets. This may indicate that there is something unique 
about the modification stage of complex crater formation on Mars, 
perhaps as a result of significant portions of both rock and water ice 
in the crust. 


It is difficult to determine exactly what is unique about Mars be- 
cause we don't really understand why a dependence on c/pg fits the 
data so well. While a dependence on hydrostatic pressure may seem 
to make sense for the initial collapse, the final crater depth, and the 
depth onset of terracing, there is no obvious reason why the trend 
applies to peak ring onset. Furthermore, why collapse occurs at all is 
somewhat puzzling, as interior slopes in simple craters are generally 
below the angle of repose (Melosh, 1977). Our results from the non- 
linear models further complicate the issue by indicating that the depen- 
dence on c/pg need not be linear; exaggerating or minimizing strength 
differences can offset a lower or higher exponential dependence on 
c/pg. However, because the trends are all roughly parallel in log-log 
space, values of b are required to be roughly similar for all the mea- 
surement types, perhaps within ±0.5. 

It is perhaps indicative of how little is really understood about 
complex crater formation that very few models in the literature make 
quantitative testable predictions of interplanetary variations in crater 
shape, and those few do not fit the data. Our inversions confirm the 
long held qualitative belief that hydrostatic pressure vs. crustal 
strength is important in complex crater formation and in fact seem 
to indicate that c/pg is the dominant factor. A quantitative physical 
model that explains in detail this apparent importance has yet to be 
developed. The general inversion techniques we present can be used 
as powerful tools for testing future models. 
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